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Summary 

A technique is presented for the calibration of a 
hemispherical- tipped, five-hole probe having a 0.125- 
in. diameter. Equations are derived from the poten- 
tial flow over a sphere relating the flow angle and 
velocity to pressure differentials measured by the 
probe. The technique for acquiring the calibration 
data and the technique for calculating the calibra- 
tion coefficients are presented. The accuracy of the 
probe in both the uniform calibration flow field and 
the nonuniform flow field over a 75°-swept delta wing 
is discussed. 


Introduction 

The measurement of flow field quantities, such 
as flow velocity and angularity, has been conducted 
with the use of pressure probes with great success for 
many years. Through the years, a variety of pressure 
probe designs and calibration techniques have been 
developed. (See ref. 1.) The hemispherical-tipped, 
five-hole probe is a good probe for the measurement 
of flow velocity and angularity in low-speed flows 
because of its high sensitivity to flow angle and its 
insensitivity to Reynolds number based on probe 
diameter R j (i.e., 2 x 10 3 < < 150 x 10 3 ). 

This report details the theoretical derivation of 
calibration equations and the calibration technique 
for a hemispherical-tipped, five-hole probe having a 
0.125-in. diameter. This probe is used in the Langley 
Basic Aerodynamics Research Tunnel (BART). The 
potential flow over a sphere is used to derive equa- 
tions relating total-flow angle, roll angle, and velocity 
to pressure differences. These equations are general- 
ized and experimental data are then fit, in a least- 
squares sense, to obtain the calibration coefficients 
used to measure the flow over the 75°-swept delta 
wing. This calibration technique is an extension of 
the method used in reference 2, and the theoretical 
derivation of the calibration equations is also based 
on reference 2. 
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fi 


calibration constants defined 
by equations (19), (20), (21), 
and (22) 

local span of 75°-swept delta 
wing, ft 

particle diameter, ft 

probe diameter, ft 

dimensionless pressure on 
port i defined as (p, — p s )/q 
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Pi 
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rms 

S a 
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Us 

Uoo 

U, V, w 
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x,y,z 
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a 0 


ft 

fto 


root chord of 75°-swept delta 
wing (1.866 ft) 

Stokes number, jg ^ Us ®p Ebt 

pressure at port i, lb/ft 2 

static pressure of flow, lb/ft 2 

difference in pressure between 
pitch ports, lb/ft 2 

difference in pressure between 
yaw ports, lb/ft 2 

dynamic pressure of flow, 

hpuL ib/ft 2 

quantity defined by equa- 
tion (15), lb/ft 2 

Reynolds number based on 
root chord, VL/v 

Reynolds number based on 
probe diameter, Vdjv 

vortex core radius, ft 
root mean square 
sensitivity of probe to pitch, 



sensitivity of probe to yaw, 

swirl velocity near vortex core, 
ft /sec 

free-stream velocity of flow, 
ft /sec 

streamwise, lateral, and 
vertical velocity components, 
respectively, in body axis 
system, ft /sec 

velocity of flow at surface of 
sphere, ft /sec 

distance from coordinate origin 
in tunnel axis system, ft 

pitch angle of total-pressure 
port to stagnation point, deg 

pitch error due to manufactur- 
ing, deg 

yaw angle of total-pressure 
port to stagnation point, deg 

yaw error due to manufactur- 
ing, deg 


e 


total angle of total-pressure 
port to stagnation point, deg 

0{ angle between the ith port and 

stagnation point, deg 

00 total angle where q p becomes 

zero, singularity point of 
equation (19), deg 

v kinematic viscosity, ft 2 /sec 

p density, slugs/ft 3 

0 roll angle at total-pressure 

port from lower a port to 
stagnation point, deg 

0 o fixed roll-angle error between 

probe calculation and sting, 
deg 

ip yaw angle of probe, deg 

00, pitch misalignment in pitch between 

probe and free stream, deg 

*00, yaw misalignment in yaw between 

probe and free stream, deg 

Subscripts: 
a air 

meas measured 

p particle 

Calibration Technique 

The probe to be calibrated is a yaw-pitch probe 
with a total-pressure port at the forward point of a 
hemispherical tip and a ring of six interconnected 
static ports approximately eight probe diameters 
from the tip. Details of the probe are presented in 
figure 1. To measure the angle of the velocity vec- 
tor, four ports are placed at approximately 45° to the 
total-pressure port in the directions of yaw and pitch. 
The total-pressure port and static-pressure port are 
numbered ports 1 and 2, respectively. (See fig. 2.) 
When the probe is lined up with the local velocity 
vector, the difference between these pressures gives 
the standard incompressible measurement of the dy- 
namic pressure. When looking upwind, the right and 
left ports are called the /? ports and are labeled 3 
and 4, respectively. These will give the angle of yaw, 
which is the angle /?. The top and bottom ports are 
called the a ports and are labeled 5 and 6, respec- 
tively. These will give the angle of pitch, which is the 
angle a. 

To specify the orientation of the velocity vector 
with respect to the probe, the usual coordinates are 


the angles a and (3 , which are rotated about the y- 
and 2r-axes of the probe as shown in figure 2. For 
a complete solution of the calibration problem to 
large angles, it is more convenient to use the angles 0 
and 0. Since the pressures on the hemispherical tip 
are assumed to be symmetrical about the stagnation 
point, the significant angle is the total angle 0 that 
the probe makes with the velocity vector. To spec- 
ify the direction of the velocity vector relative to the 
probe, a polar angle 0 is taken about the probe and 
referenced to the lower a port. Spherical trigonom- 
etry yields the following conversion between the two 


coordinate systems: 

tan a = tan 0 cos 0 ( la) 

sin /? = sin 0 sin 0 (lb) 

cos 0 = cos a cos f3 (lc) 

tan 0 = tan /?/ sin a (Id) 


Calibration Based on the Potential Flow 
Over a Sphere 

This calibration technique closely follows the 
method described in reference 2. The pressure on a 
sphere is a function only of the total angle 0 from the 
stagnation point. A nondimensional quantity f{ may 
be defined as the pressure at port i (p t ) minus the 
static pressure (p s ), divided by the dynamic pressure 
(»): 

At*.') = ~ (2) 

This quantity can be determined theoretically by 
using the potential-flow theory for uniform flow over 
a sphere. If axisymmetric flow is assumed, the 
superposition of a free stream and a doublet flow 
gives the velocity on the surface of the sphere as 

V{0) = \u o osin0 (3) 

Bernoulli’s equation is then used to obtain the pres- 
sure distribution on the sphere as 

Pi + \p y2 =Ps + \pulc (4) 

By using this equation and the definition of /^, one 
obtains 

/i(0»-) = J cos 2 6i - ^ (5) 

The most important quantities are the differences 
in pressure between the two a ports and the two /? 
ports. These are labeled A p a and A respectively, 
with the convention that they are positive in the 
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a- and /3-directions. From the definition of / t -, the 
following relations are obtained: 

A p a = q{fe ~ h) (6) 

Ap/3 = q{h - /*) (7) 

By using the functional form of /;, one obtains the 
relations 

Ap a = 7 q (cos 2 $6 - cos 2 05) (8) 

4 

App = j<?(cos 2 03 - cos 2 04) (9) 

These relations can be converted to functions of 0 
and (f> by reference to figure 2. The law of cosines for 


spherical triangles yields 

cos 0g = cos 0 cos 45° -h sin 0 sin 45° cos </> ( 10a) 

cos 05 = cos 0 cos 45° — sin 0 sin 45° cos </> (10b) 

cos 04 = cos 0 cos 45° — sin 0 sin 45° cos <j> (10c) 

cos 03 = cos 0 cos 45° + sin 0 sin 45° cos <f> ( lOd) 


Substituting equations (5) and (10) into (15) 
reduces equation (15) to 

Qp = q (% cos26+ ^ ( 16 ) 

Substituting equation (16) into equation (14) 
yields 

(Ap 2 + Apjg) 1 / 2 | sin 20 

q P || cos 20 + Jtj 

It should be noted that this has a singularity at 
54.7°. The technique given here does not use the 
measurements made on the static-pressure ports, 
which can have significant errors at large angles to 
the velocity. However, once the dynamic pressure q 
and the total angle 0 are known, the static pressure 
can be calculated from the quantity f\ by using the 
following formula: 

Pi -Ps = g(^cos 2 0 - j) (18) 

Equations (13), (16), (17), and (18) determine the 
primary quantities from which the angles 0 and <f> and 
the pressures q and p s can be found. 


By substituting these equations into the equations 
for A p a (eq. (8)) and App (eq. (9)), one obtains 

9 

A p a = —q sin 20 cos (f> (11) 
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App = -q sin 20 sin <f> (12) 

The dependence on 0 and </> can be separated by 
taking the ratio and the square root of the sum of 
the squares of these two equations: 

= tan <f> (13) 

APa 

(Ap 2 + Apjg) 1 / 2 = ^<7 sin 20 (14) 

The square-root quantity may be nondimension- 
alized by dividing by the quantity q p , defined as the 
pressure on port 1 minus the average of the pressures 
on the four angle ports. This relation can be written 
in terms of as 


Generalization of Calibration Technique 

For a number of reasons, the potential-flow cali- 
bration may not be satisfactory for a given yaw-pitch 
probe. The decrease in pressures with 0 given by 
the potential-flow theory is ideal for a sphere only 
while the probe is a hemispherical-tipped cylinder, 
and thus the decreased pressures differ from experi- 
mental values at large distances from the stagnation 
point. The placement of the ports may be in error 
because of manufacturing problems. This means that 
there will be fixed errors in the a - and /3-directions, in 
the rotation angle </>, and in the sums and differences 
of the pressures on the a and /3 ports. For these 
reasons the theoretical calibration is generalized to 
include some experimentally determined parameters. 
While maintaining the same form, the constants of 
the 0 dependence are made arbitrary so that equa- 
tions (16), (17), and (18) are generalized to become 

(Ap 2 + Apg) 1 / 2 ^ 8in 20 

q p cos 20 — cos 20o 

q p = q [C(cos 26 — cos 26q) + D( cos 6 — cos 0q)] (20) 


q P = q 


h - ^(/3 + h + h + h) 


(15) 


pi - p s = q[E cos 2 0 + FcosO + G) (21) 


3 


The constant 0 q 1S the singularity point of equa- 
tion (19) and must be the same for equation (20) 
to determine correct values for q. The calculation 
for <j) (eq. (13)) is made arbitrary by the subtraction 
of a constant (j ) o and the multiplication by a constant 
B as follows: 

= tan(<£ - <f> 0 ) (22) 

Ap a 

It should be noted that the constants B and do 
not appear in the generalization of reference 2. 

Determination of Calibration Constants 

The experimental technique for acquiring the five- 
hole probe calibration data and the calculation of the 
calibration constants are described in this section. 
The 0.125-in-diameter five-hole probe was mounted 
in the test section of the Langley Basic Aerodynamics 
Research Tunnel on a C-strut so that the probe tip 
was always positioned at the same location when 
the probe was yawed. Pressure data were obtained 
using an electronic scanning pressure system with 
±1 lb/in^ transducers. Four sets of calibration data 
were obtained at eight free-stream velocities. A set 
of calibration data was obtained by yawing the probe 
between -90° and 90° in 1° increments. The four 
data sets were obtained by rolling the probe through 
90° increments. (Positive roll is counterclockwise 
when looking downstream at the probe tip.) A 
typical data set is presented in figure 3 for a free- 
stream velocity of 160 ft/sec and a roll of 0°. These 
data sets were then used to calculate the calibration 
coefficients. The procedure used to calculate the 
calibration coefficients involved six computer codes 
and is described below. 

Before the discussion of the methods used to 
calculate the probe errors and calibration coefficients, 
a few comments about the experimental procedure 
are necessary. When the probe is rolled 0° or 180°, 
the data acquired are for the calibration of the f3 
ports (ports 3 and 4). When the probe is rolled 90° or 
270°, the data acquired are for the calibration of the 
a ports (ports 5 and 6). The probe is always yawed 
in the same direction; therefore, when the probe is 
rolled 0°, 

P = +0 

when rolled 90°, 

a = —ifi 

when rolled 180°, 

P = -0 

and when rolled 270°, 

a = +\j) 


There are two sets of values for the constants 
0 q, 00? A, B,C, D,E,F, and G. One set is derived 
for the a ports and the other set for the f3 ports. 
Each set of constants is used to calculate a value of 
6 and 0 from the pressures measured by the five-hole 
probe. These two independent measures of 6 and 0 
are converted into a and /3 and then are combined 
by using a cosine weighting function based on the 
roll angle 0. The final values of 6 and (f) are obtained 
using equations (lc) and (Id), respectively. 

Step 1 — Calculation of probe manufacturing errors 
and misalignment errors between the probe and the free 
stream . The first step in the probe calibration proce- 
dure is the determination of the probe manufactur- 
ing errors and the misalignment errors between the 
probe and the free stream. These errors are defined 
as <*o and /3q (the errors caused by improper place- 
ment of the pressure orifices) and 0o,yaw and 0o, pitch 
(the misalignment errors of the probe with the free 
stream). Roll data that are 180° out of phase are 
used to calculate these errors. 

To understand how the misalignment error of 
the probe with the free stream in yaw (0o,yaw) IS 
calculated, assume that the probe is made perfectly 
(i.e., no probe error) and that we are looking at data 
for the (3 ports (rolls of 0° and 180°). If the probe is 
misaligned in yaw by some positive increment, then 
the probe reading A pp at zero yaw will be negative 
for the roll of 0° and positive for the roll of 180° 
(fig. 4(a)). The probe will have the same magnitude 
of sensitivity to flow angle in either roll; however, 
the sign of the sensitivity will be the opposite for 
the roll of 180°. Plotting the two roll data sets 
on the same plot yields figure 4(b), and the probe 
misalignment error is found at the intersection of the 
two curves. One should note that this error should 
be the same for both the f3 calibration data sets and 
the a calibration data sets because of the method of 
obtaining the calibration data. 

Probe error due to manufacturing is found by 
a similar method. Assume that the probe is per- 
fectly aligned with the free-stream direction, but the 
(3 ports have been incorrectly placed as illustrated 
in figure 5(a). For the roll of 0°, the five-hole probe 
would measure some flow angle at 0° yaw (+Ap^). 
For the roll of 180°, the five- hole probe would mea- 
sure the same value of +A pp. Once again the sensi- 
tivity of the probe to flow angle would have the same 
magnitude but would differ in sign for the two roll 
cases. If one uses this information, a figure can be 
constructed that allows the probe manufacturing er- 
ror to be calculated (fig. 5(b)). This technique is also 
used for the a ports. 
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Because the probe is not rotated in pitch, the 
only method for calculating the misalignment of the 
probe in pitch is by using the flow angle predicted by 
the probe and the probe error due to manufacturing. 
Figure 6 presents sketches that help in the under- 
standing of this pitch misalignment. The measured 
A p a is a function of the error in manufacturing of 
the a ports and the misalignment of the probe with 
the free stream in pitch so that 

V’O, pitch = [(&Pa/q)/S a ] - ao 

for the rolls of 0° and 180°, where S a is determined 
by the slope of the calibration data. (Determination 
of Sp is demonstrated in fig. 5.) 

Figure 7 presents the runs for rolls of 0° and 180° 
for the five-hole probe with a free-stream velocity of 
160 ft/sec and shows the combined effect of probe er- 
ror due to manufacturing and misalignment. These 
errors can be measured directly off the plot; however, 
the code used for step 1 of the probe calibration cal- 
culates the errors by fitting a least-squares linear fit 
to the data in the yaw range from —12° to 12° and 
calculates the intersection of the two curves with each 
other and with the horizontal axis. These intersec- 
tions yield the information required to deduce the 
probe manufacturing and misalignment errors. One 
last point, because the probe is rolled, the combina- 
tion of probe errors due to manufacturing and mis- 
alignment with free stream is different for each roll 
angle to correct the calibration data properly. Thus, 
for a roll of 0° , 

0 — V 7 ~ V’O, yaw — 00 
a ~ “V’O, pitch “ <*0 

for a roll of 90°, 

0 — “V’O, pitch “ 00 
a = -(l/> - V’O, yaw + ao) 

for a roll of 180°, 

0 = “(V’ “ V’O, yaw + 0o) 

<x = “V’O, pitch + a 0 

and for a roll of 270°, 

0 = “V’O, pitch + 00 

Ot — Tp V’O, yaw “ 


Step 2 — Calculation of calibration coefficients 0q 
and A . Once the probe manufacturing errors and the 
probe misalignment with the free stream are sub- 
tracted from the calibration data, the calibration co- 
efficients 0 q and A are found by fitting the calibra- 
tion data to equation (19). Although yaw data are 
obtained between —90° and 90°, only yaw data be- 
tween —45° and 45° are used for the curve fit. The 
data are fit in a least-squares sense by using the dif- 
ferential correction method detailed in reference 3. 
It is assumed for the first iteration that B and <f>o 
are equal to their theoretical values (1 and 0, respec- 
tively). 

Step 3 — Calculation of calibration coefficients <f o 

and B. Once 6q and A are obtained from step 2, 
the calibration coefficients </>o and B are found by 
fitting the calibration data to equation (22) in a least- 
squares sense by using the differential correction 
method. Once these constants have been found, they 
are used in the step 2 calculations to yield a better 
fit for 6q and A. Steps 2 and 3 are iterated until the 
constants A, E, and 4>q have converged (usually 
requiring three iterations). The iteration technique 
is required by the software used to calculate the yaw 
and pitch from the differential pressures measured by 
the five-hole probe. (See the appendix.) 

Step 4 — Calculation of calibration coefficients C 

and D . Once the constants A, B , and V>o have 
converged, they are used to calculate the calibration 
coefficients C and D. These coefficients relate the 
pressures measured at the tip of the probe to the 
local dynamic pressure of the flow and are found in 
a manner similar to the others. The calibration data 
are fit in a least-squares sense to equation (20) by 
using the differential correction method. 

Step 5 — Calculation of calibration coefficients E, 

F 9 and G . The last step in the probe calibration is the 
calculation of the calibration coefficients F, F, and G. 
These coefficients relate the pressures measured on 
the probe tip to the static pressure in the flow. The 
best-fit values for the calibration coefficients V’O) 
A , B y G, and D are used while the calibration data 
are fit, in a least-squares sense, to equation (21) by 
using the differential correction technique mentioned 
previously. 

Technique for Measuring Unknown Flow 
Velocity and Angularity 

Once the calibration coefficients are obtained, 
their dependence on free-stream velocity (Reynolds 
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number) can be assessed. Figure 8 presents the cali- 
bration coefficients A, B, C , and D as a func- 

tion of free-stream velocity. (This set of calibration 
data did not allow for the calibration of the static 
pressure.) A third-order polynomial curve fit was 
made to each of the calibration coefficients, and these 
polynomials were incorporated into the subroutine 
that calculates the flow angularity and velocity from 
the measured pressures. This subroutine is presented 
in the appendix. The flow angularity and velocity 
were found in an iterative manner. The free-stream 
velocity is used to calculate the first set of calibra- 
tion coefficients, and on the subsequent iterations the 
velocity predicted by the five-hole probe is used to 
calculate the calibration coefficients. The flow angu- 
larity and velocity are found when the velocity pre- 
dicted by the probe has converged to one value. 

Accuracy of Five-Hole Probe 

The measured flow angle and velocity have been 
compared with the known values of flow angle and 
velocity for the calibration data. The probe error 
was calculated by taking the root mean square (rms) 
of the difference between the known and the mea- 
sured flow angle and velocity. These rms errors are 
presented as the solid lines in figure 9. The figures 
indicate that the probe measurements for 6 < 50° 
are as follows: the pitch angle a within 0.5°, the 
yaw angle /3 within 1.5°, and the dynamic pressure q 
within 0.02 lb/ft 2 . Figure 9 also presents the effect 
of Reynolds number for the range of Reynolds num- 
ber that the probe was calibrated over. These data 
(the dashed lines in fig. 9) were found by taking the 
calibration constants for the highest Reynolds num- 
ber and applying them to the data acquired for the 
lowest Reynolds number and vice versa. These data 
indicate little effect due to Reynolds number in the 
measurements of a and /3; however, the error in the 
measurement of the dynamic pressure increases by a 
factor of 3. 

These characteristics were measured with a uni- 
form onset flow. When the probe is introduced into 
the flow over a model, the flow field seen by the 
probe is no longer uniform. The gradients in the 
flow field can affect the accuracy of the probe and 
the probe may disturb the flow field. To help under- 
stand the magnitude of the errors caused by the five- 
hole probe in a nonuniform flow field, a comparison 
has been conducted (ref. 4) between flow field data 
obtained over a 75°-swept delta wing with the five- 
hole probe and data obtained by a three-component 
laser Doppler velocimeter (LDV). The LDV is ca- 
pable of obtaining accurate velocity measurements 
in flow fields with reverse flows, large shear gradi- 
ents, and velocity fluctuations. The BART LDV 


system is a three-color, orthogonal, cross-fringe con- 
figuration with the receiving optics package mounted 
90° off-axis. The 514. 5-, 496.5-, and 476.5-nm wave- 
lengths are used to measure the lateral (u), stream- 
wise (u), and vertical (w) velocity components, re- 
spectively. Bragg cells are used to provide directional 
measurement capability in all three velocity compo- 
nents. The sample volume is spherical in shape and 
has been calculated to be approximately 150 fim in 
diameter. The optics and laser move as a unit on a 
traversing system that provides 1 m of travel, with 
10 /i m of resolution, in all three axes. 

The flow field was seeded with 0.8-pm polystyrene 
latex microspheres. The seed particles were sus- 
pended in a mixture of alcohol and water and were 
injected into the flow upstream of the honeycomb by 
using an atomizing spray nozzle. Typically, 500 to 
4096 velocity samples were obtained at each measure- 
ment location in the flow field. The actual number of 
samples depended on the particular location in the 
flow field and the particle seeding rate. 

The ability of a particle to track the streamlines 
in the flow field, and thus the accuracy of the LDV, is 
directly related to the size of the particle. Theoretical 
predictions of particle trajectories in various flows 
were reported in references 5 through 7. Dring and 
Suo (ref. 5) concluded that the particle trajectory in a 
free- vortex swirling flow is governed primarily by the 
Stokes number (7Vg t ); and when the Stokes number 
is less than 0.01, the particle will follow the circular 
streamlines of the free vortex. 

The O.S-fim particles used during this test have a 
density p p of 2.03727 slugs/ft 3 . The Stokes number 
for the particles is 0.007, based on the radius and the 
swirl velocity at the edge of the vortex core. The 
numerical procedure described by Dring (ref. 7) was 
used to predict the particle trajectories for the vor- 
tices that were measured during this investigation. 
The predictions show that the particles used during 
this test will follow the streamlines of the vortex with 
an accuracy of about 1 percent. 

The velocity surveys were made over a 75°-swept 
delta wing at an angle of attack of 20.5° and a 
Reynolds number of 1 X 10 6 based on the root chord. 
These velocity surveys are presented in figure 10. 
The differences in the u, u, and w components of the 
flow field are presented in figure 11. In the core of the 
vortex the five-hole probe has errors in the u, u, and 
w components of velocity of 25, 17, and 35 percent, 
respectively. The errors were calculated assuming 
that the LDV measurements were the reference. The 
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equation used to calculate the u-component error is 


terror — 100 

( U LDV 


U LDV ~ U FH 

+ U LDV + ^LDv) ^ 


where the subscript FH denotes the five-hole probe. 
The errors in v and w are calculated in a similar 
fashion. 

Figure 12 presents the gradients in the v and 
w components of velocity measured by the five-hole 
probe for the same data set. A comparison of fig- 
ures 11 and 12 shows that the probe does a reason- 
able job of measuring the flow field quantities in re- 
gions of low gradients (with a probe error less than 
5 percent when the velocity gradient is less than 
800^ /ft). However, in the vortex core, very high 
gradients (> 3000^/ft) yield differences of more 
than 20 percent between the two measurement tech- 
niques. 

Concluding Remarks 

A technique has been presented for the calibration 
of a hemispherical-tipped, five-hole probe with a 
0.125-in. diameter. Equations were derived from 
the potential flow over a sphere relating the flow 
angle and velocity to pressure differentials measured 
by the probe. Four sets of calibration data were 
acquired at eight free-stream velocities. A set of 
calibration data was obtained by yawing the probe 
between —90° and 90° in 1° increments. The four 


data sets were obtained by rolling the probe through 
90° increments. The calibration data are fit to the 
equations derived from the potential flow over a 
sphere to obtain the calibration coefficients, which 
are then used to convert measured pressures into flow 
angularity and velocity data. 

To assess the ability of the five-hole probe to mea- 
sure flow angularity and velocity, the known calibra- 
tion data were compared with those measured by the 
five-hole probe. For the uniform flow field measured 
for the probe calibration, the probe measurements 
were as follows for a total angle of total-pressure 
port to stagnation point less than 50°: the pitch an- 
gle within 0.5°, the yaw angle within 1.5°, and the 
dynamic pressure within 0.02 lb/ft 2 . For the non- 
uniform flow over a 75°-swept delta wing, com- 
parisons between the five-hole probe and a three- 
component laser Doppler velocimeter showed that 
the probe performed a reasonable job of measuring 
the flow field quantities in regions of low gradients 
(with a probe error less than 5 percent when the ve- 
locity gradient is less than 800^ /ft). However, in 

the vortex core, very high gradients (> 3000^/ft) 
yielded differences of more than 20 percent between 
the two measurement techniques. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
September 20, 1988 


Appendix 


Listing of Five-Hole-Probe Subroutine 

This subroutine is used to convert the pressures measured by the yaw-pitch probe into flow velocity and 
angularity. 

SUBROUTINE FIVE_HOLE(P, ALPHA, BETA, THETA, PHI, Q, PS, VEL.PREV) 

C 

C FIVE-HOLE CONVERTS MEASURED PRESSURES INTO FLOW ANGLES 
C AND VELOCITY DATA USING THE CALIBRATION EQUATIONS DERIVED 
C IN APPENDIX A OF REFERENCE 2 

C SUBROUTINE WRITTEN BY KJELGAARD IN AUGUST 1986 

C DESCRIPTION OF VARIABLES 
C INPUT 

C P - ARRAY CONTAINING PRESSURES MEASURED BY 5-HOLE PROBE 
C P(l) CONTAINS PRESSURE AT PORT 1 

C P(2) CONTAINS PRESSURE AT PORT 2, ETC. 

C OUTPUT 

C ALPHA - PITCH ANGLE MEASURED BY PROBE 
C BETA - YAW ANGLE MEASURED BY PROBE 

C THETA - TOTAL FLOW ANGLE MEASURED FROM PORT 1 TO STAGNATION 
C POINT 

C PHI - FLOW ROLL ANGLE MEASURED FROM LOWER ALPHA PORT TO 
C STAGNATION POINT 

C Q - DYNAMIC PRESSURE OF THE FLOW AT THE PROBE TIP 

C PS - STATIC PRESSURE MEASURED AT THE PROBE TIP 

C VEL.PREV - VELOCITY USED TO CALCULATE PROBE CONSTANTS 
C 

REAL P(*) 

REAL ALPH (2) ,BET(2) ,THET(2) ,QT(2) ,PST(2) 

REAL A(2),THETA0(2),PHI0(2),B(2),C(2),D(2),E(2),F(2),G(2) 

EXTERNAL POLY 
DPA=P(6)-P(5) 

DPB=P(3)-P(4) 

QP=P(1)-.25*(P(3)+P(4)+P(5)+P(6)) 

A1=SQRT(DPA**2+DPB**2)/QP 

C 

C PROBE CONSTANTS 

C 

A(1)=POLY(VEL_PREV,.13930,-2.5727E-3,1.43414E-5,0) 

A(2)=POLY(VEL_PREV, 3.1607, -8.1927E-3,7.1750E-5,-1.8642E-7) 

THETAO(1)=POLY(VEL_PREV, 55.595, -5.8621E-2,4.7727E-4,-1.1427E-6) 
THETAO(2)=POLY(VEL.PREV, 57.134, -9.3249E-2,8.0274E-4,-2.033lE-6) 

B(l)=POLY(VEL_PREV,. 31547, 2.4019E-2,-2.0699E-4,5.5156E-7) 

B(2)=POLY(VEL_PREV,. 67481, 9.7006E-3,-7.1074E-5,1.4787E-7) 
PHI0(l)=POLY(VEL_PREV,-.48446,-4.9767E-2,4.5553E-4,-1.2265E-6) 

PHIO(2)=POLY(VEL_PREV, -4.2108, 4.5463E-2,-2.7104E-4,4.4549E-7) 
C(l)=POLY(VEL_PREV,.25419,-1.6741E-3,9.3270E-6,-8.6405E-9) 

C(2)=POLY(VEL_PREV,. 38568, -5.2703E-3,3.7782E-5,-8.2944E-8) 
D(l)=POLY(VEL_PREV,.74529,1.7064E-2,-1.1424E-4,2.2726E-7) 

D(2)=POLY(VEL_PREV, . 23948, 2.9862E-2,-2.1885E-4,5.0793E-7) 

C 

C ITERATIVELY SOLVE FOR THETA 


c 

DO 1 1=1,2 
THET(I)=0. 

THSTEP=10. 

C 

C CALCULATE PHI 
C 

IF (DPA .EQ. 0) DPA=.00001 
PHIM=ATAND(B(I)*DPB/DPA) 

IF (DPA .LT. 0) PHIM= 180.+PHIM 
IF (PHIM .LT. 0) PHIM=PHIM+360. 

PHIT=PHIM-PHIO(I) 

A2=l 

DO 2000 ILOOP=1,50 
IF (ILOOP .EQ.50) THEN 

* WRITE(6,‘(“ TOO MANY ITERATIONS IN CALB”)’) 

THET(I)=.5*ACOSD(A(I)*SIND(2.*THETAO(I))/A1+ 

* COSD(2.*THETA0(I))) 

GO TO 3170 

END IF 

THET(I)=THET(I) + THSTEP 
A3=COSD(2.*THET(I))-COSD(2.*THETAO(I)) 

IF (A3.EQ.0.0) THEN 
WRITE(6,‘(“ PROBLEM IN CALIB”)’) 

WRITE(6,‘(4F 10.5) ’) THET(I) , THSTEP 
A3=lE-5 
END IF 

A2=A1-A(I)*SIND(2.*THET(I))/A3 
IF (A2 .LT. 0) THEN 
THET(I)=THET(I)-THSTEP 
THSTEP=THSTEP*.5 
GO TO 3160 
END IF 

IF (A2 .LE. .02) GO TO 3170 
3160 CONTINUE 
2000 CONTINUE 

3170 ALPH(I)=ATAND(TAND(THET(I))*COSD(PHIT)) 
BET(I)=ASIND(SIND(THET(I))*SIND(PHIT)) 
QT(I)=QP/(C(I)*(COSD(2.*THET(I))-COSD(2.*THETAO(I)))+ 

* D(I)*(COSD(THET(I))-COSD(THETAO(I)))) 
PST(I)=P(l)-QT(I)*(E(I)*COSD(2.*THET(I))+F(I)*COSD(THET(I)) 

* +G(I)) 

1 CONTINUE 

C 

C USE PHI FOR WEIGHTING FUNCTION FOR COMBINATION OF ALPHA 
C AND BETA CONSTANTS AND COMBINE 
C 

WGHT=COSD (PHI) * *2 
WGHT1=1-WGHT 

BETA=BET(1)*WGHT+BET(2)*WGHT1 
ALPHA=ALPH ( 1) * WGHT+ALPH (2) * WGHT1 
Q=QT(1)*WGHT+QT(2)*WGHT1 
PS=PST(1)*WGHT+PST(2)*WGHT1 
C 

C USE THESE TO CALCULATE THETA AND PHI 



THETA=ACOSD(COSD(ALPHA)*COSD(BETA)) 
IF (ALPHA .EQ. 0) THEN 
PHI=SIGN(90.,BETA) 

ELSE 

PHI=ATAND(TAND(BETA)/SIND(ALPHA)) 
END IF 

IF (ALPHA .LT. 0) PHI= PHI+180. 

IF (PHI .LT. 0) PHI= PHI+360. 

RETURN 

END 

FUNCTION POLY(V,A,B,C,D) 

POLY = A + B*V + C*V*V + D*V*V*V 

RETURN 

END 
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Figure 1. Details of five-hole probe. All linear dimensions are in inches. 










Figure 8. Variation of calibration coefficients with free-stream velocity. 






.11s of 0° and 180°. 



Figure 8. Concluded. 
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(a) Crossflow velocity vectors measured with five-hole probe. 



(b) Crossflow velocity vectors measured with three-component laser Doppler velocimeter. 

Figure 10. Velocity surveys obtained over 75°-swept delta wing at an angle of attack of 20.5 , i/L — 0.9, and 
R = 1 x 10 6 . 
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